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ABSTRACT 


In this paper the crack problem of a shallow shell with two nonzero 
curvatures is considered. It is assumed that the crack lies in one of 
the principal planes of curvature and the shell is under Mode I loading 
condition. The material is assumed to be specially orthotropic. After 
giving the general formulation of the problem the asymptotic behavior of 
the stress state around the crack tip is examined. The analysis is based 
on Reissner's transverse shear theory. Thus, as in the bending of cracked 
plates, the asymptotic results are shown to be consistent with that ob- 
tained from the plane elasticity solution of crack problems. Rather 
extensive numerical results are obtained which show the effect of material 
orthotropy on the stress intensity factors in cylindrical and spherical 
shells and in shells with double curvature. Other results include the 
stress intensity factors in isotropic toroidal shells with positive or 
negative curvature ratio, the distribution of the membrane stress resultant 
outside the crack, and the influence of the material orthotropy on the 
angular distribution of the stresses around the crack tip. 

1. Introduction 

The crack problem in shallow shells by using a transverse shear 
theory [1,2] has previously been considered for cylindrical and spherical 
shells only [3-5]. The results given in [3-5] as well as those obtained 
from the classical shell theory (e.g., [6-8]) indicate that due to the 
curvature effects the stress intensity factors in shells may be consider- 
ably higher than that in flat plates having the same crack length and the 
same thickness. Also, the results given in [5] and [8] show that, unlike 
the infinite plate problem, the material orthotropy may have a significant 
influence on the stress intensity factors in shells. Therefore, from the 
viewpoint of practical applications it does seem to be important to study 
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the problem for shell geometries other than cylinders and spheres and 
also to consider the effect of material orthotropy. 

In this paper the basic problem of a shallow toroidal shell with 
two unequal and nonzero curvatures is considered. It is assumed that the 
material is specially orthotropic, the through crack is located in one 
of the principal planes of curvature, and the plane of the crack is a 
plane of symmetry with respect to loading as well as to the geometry of 
the shell. To remove the inconsistency in the asymptotic behavior of the 
bending and transverse shear resultants, a transverse shear theory is used 
in the analysis [1,2] (see the discussion given in [3-5] for shells and, 
for example, [9-11] for plates). 


2. The Integral Equations 

The geometry of the shallow shell under consideration is shown in 
Fig. 1. If the material is orthotropic (with x-| and x 2 as the axes of 
orthotropy), the engineering material constants are defined through the 
following stress strain relations: 

e ll = r i (a ll" v l 0 22 ) * e 22 = r 2 ^ a 22" v 2 a ll^ » 

e-j 2 - 0 i 2^®12 * v l^l - v 2^2 * ) 


With the four independent material constants the differential operators 
arising from the formulation of the shallow shells do not seem to be 
factorable and consequently, the analysis becomes intractable. However, 
it can be shown that if the elastic constants are related through 





2(1 +1^ v 2 ) 


( 2 ) 


then, with a simple coordinate transformation, the orthotropic shell 
equations may be reduced to essentially that of isotropic shells (see, 
for example, [8] and [3-5]). The factorization condition (2) implies that 
the shell has only three independent elastic constants. Such a material 
(in plate or shell form) is said to be specially orthotropic. If we now 
define 
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E = /EjE^ 


= , G ay = E/2(l+v) , c = (E^Eg) 


the stress-strain relations (1) become 


G 11 ^'^■■ vo 22 ) » e 22 = l^ c ^ a 22 ” va l 1 ^ * e 12 = 25“, 


where E and v are the effective modulus and Poisson's ratio, c is the 
stiffness ratio, and from (2) and (3) it follows that in the specially 
orthotropic materials the measured shear modulus G^ is equal to the 
(calculated) effective shear modulus G_„ . 

a V 

Referring to Fig. 1 and the normalized quantities defined in Appendix 
A, and to [3-5] for details, in terms of a stress function <(> and the (out 
of plane) displacement component w the shell problem may be formulated 
as follows: 

v% - v x w = 0 » ( 5 ) 

v 4 w + x 2 (l-<v 2 )v 2 <j> = 0 , (6) 

o 

<Vi{i-tp-w = 0, (7) 


where 


v 2 n - a = 0 , 


2 _ 2 d ? r . ,2 a 

x ~ i ^ 7 ” 12 3x3 y 2 ^7 * 






The shell parameters Xj, x^» and k are defined in Appendix A, 

6 and 8 are the angles of rotation of the normal to the shell surface, 
x y 

and the curvatures of the shell are defined by 


1 _ a 2 Z 1 a 2 Z 1 _ a 2 Z 


ax 2 * «2 


ax| * R 12 


ax-jaxg * 


where Z = Z(x^,X2) is the equation of the middle surface of the shell 
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x 1 and x 2 being the coordinates in the tangent plane. The normalized 
stress, moment, and transverse shear resultants are given by 


2 2 2 
,i _ 3 <f> ,i . 3 ♦ m _ 3 $ 

N xx - ’ N yy - • “xy - Sw ’ 


M. 


xx 


M. 


_ __a 

hx' 

_ a 


5\3x 3y ) yy 8x 8y ) 


l-v/ 36 x + fy'j 
TT\»y *n) ’ 




V =|w + s , V -IS+B , 

x 3x x * y 3y y * 


where 


o - 3£ + 1 — v 3ft ft _ 3^ 1 — v 3ft 

p x • 3X K 2 3y * B y " 3y " K 2 3x * 


( 12 ) 

(13) 

(14) 

(15) 


The system of differential equations (5)-(8) may be solved by 
using the standard Fourier transforms. In this paper it is assumed 
that through a linear superposition the regular part of the solution 
has been separated and the problem is reduced to a perturbation problem 
in which the self-equilibrating crack surface tractions are the only 
nonzero external loads. Because of the assumed symmetry with respect 
to x 2 x 3 plane in loading and shell geometry, the stress and moment 
resultants must satisfy the following symmetry conditions: 


N xx (x,y) = N xx (-x,y), N xy (x,y) =-N xy (-x,y), 

N yy (x,y) = N yy (-x,y) , V x (x,y) =-V x (-x,y), V y (x,y) = V y (-x,y), 

M xx (x,y) = M xx (-x,y), M yy (x,y) = M yy (-x,y), M xy (x,y) =-M /y (-x,y), (16) 


It is, therefore, sufficient to consider one half of the shell only. 
Thus, in addition to the regularity conditions at infinity, the problem 
must be solved under the following boundary conditions: 


N xy (0,y) = 0, M xy (0,y) = 0, V x (0,y) = 0, -»<y<», (17) 
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08 ) 


lim N (x,y) = F,(y), |y|<v£, u(+0,y) = 0, ^c<|y|<», 

x-ho xx 1 

lim M (x,y) = F«(y), |y|<»^, B x (+0,y) = 0, v£<|y|<«. (19) 

x-*+0 xx c x 


We note that the problem (for the half shell) is one of tenth 
order. By taking Fourier transforms in y it would give ten "integration 
constants" which are functions of the transform variable and are unknown. 
Five of these unknowns must be zero because of the regularity conditions 
at x=®, three may be eliminated by using the homogeneous conditions (17) 
at x=0, and the two mixed boundary conditions (18) and (19) would give 
a pair of integral equations to determine the remaining two. We also 
note that the integral equations of this problem would be identical to 
those obtained in [5] for the spherical shell. In [5], even though 


R-j =[* 2 * because of the assumption of special orthotropy, the stiffness 
ratio c^E-j/Eg) 14 is not unity and consequently X-j fx 2 , which is also the 
case in the problem considered in this paper. Thus, by defining 


Jy u(+0,y) = G-j (y) , ^ 3 x (+0,y) = G 2 (y) , (20) 

from (18) and (19) it is seen that 


Gj(y) = 0, |y|>r€, (j-1,2) 

(21) 

fJc 

j ^_8j(y)dy = o, (j=i,2) . 

(22) 


Referring now to [5] for details the integral equations of the problem 
may be expressed as 

rfi G, (t) ri/c 2 

j -^y- dt + j l kjjty.tjGjOOdt ■ 2 1 rF 1 (y), \y\<Jc, 

" rC “vC 1 

-i 2 r/E 6p(t) ri/c 2 . _ 

•^J -f^-dt + J l k 2j (y»t)G j (t)dt = 2* J F 2 (y), \y\<Jc, (23a, b) 

X "/c "/c 1 


where the kernels k. . are known bounded functions and are given in [5]. 

' J 
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3.- The Stress State Around Crack Tips 

In the system of integral equations the interval (-/c,^ may be 
normalized by introducing the following change in variables 

t=t/v£, r)=y/Jc, 5 = x//c, s j(t) = g j ( t ) , (j-1,2) . (23) 

The integral equations may thus be expressed as 

J . p?* K ij( n ’ T )]9j( T ) dT ■ M”)* ( 1=1 > 2 ). -l<n<l . (24) 

The index of the singular integral equations (24) is +1 and their 
solution is, therefore, of the following form 

^•(x) 

g,(x) = - j— , (i=1 ,2) , ”1 <t<1 , (25) 

/l-x 2 

where h^ and h 2 are bounded functions. The system of integral equations 
(24) subject to the additional single-valuedness conditions (see (22)) 

j^g^xjdx = 0, (i=l ,2) (26) 

may be solved numerically by using the related Gaussian quadrature 
formulas [12]. 

After solving the integral equations (24) all the field quantities 
in the shell may be expressed in terms of the density functions g-j and 
g 2 as finite integrals [5]. In particular, the behavior of the stress 
state around the crack tips may be obtained by examining the corresponding 
integrals asymptotically. Referring to Appendix A, the normalized membrane 
and bending stresses are given by 


m ^ i b 1 2az ^ _m j_ _b f ^ ^ \ i , ,\ 

°ij “ N ij» °ij h~ M ij* °ij " a ij °ij* “ ^ x,y ^ * 


(27) 


If we now define the polar coordinates (r,e) at the crack tip (5=0, n=l) 
by 

5 = r sine, n-1 = r cos? , (28) 
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the asymptotic stress state at the crack tip may be expressed as 
(see [5] for details) 


a 


xx 



h 1 (l)+zh 2 (l) 

h^p+zhgO) 

2 / 2 r 

h^p+zhgO) 

2 v£r 


/5 * 0 1 _ _ 5q \ 

($• cos cos -£-) , 

/ 3 _ e.l 50\ 

cos 2 + 4 cos g ) » 

/ 1 0.1 • 50 \ 

(“5" Sln 2’ + 4' sln ”2") » 


(29a-c) 


V x = [-§■ h 1 (l)x 2 x| + h 2 (l)]i^r7? sin 0 cos|- . (30) 


Note that the asymptotic behavior of the in-plane stresses given by 
(29) is identical to that obtained from the plane elasticity solution 
of a symmetric crack problem. The transverse shear resultant given 
by (30) represents Mode III cleavage stress which is nonsingular for 
the symmetric loading under consideration. 

Referring to Fig. 1, in this problem the Mode I stress intensity 
factor at the crack tip x 2 =a is defined as follows: 


ki (X 3 ) - 1 im ✓' 2 (x 2 -a) ®'j|(O»X 2 »X 0 ) 


(31) 


Xg+a 


From the definitions given by (23), (28) and Appendix A it can 
be shown that 

F 1 X ] 

tane = -V = 7 - . 

n - 1 c Xg-a 


(32) 


If we define the polar coordinates in the x^x 2 plane (at x^=0, x 2 =a) by 


p since = x-j , p cosa = x 0 - a , 


(33) 


from (32) it follows that 


tana = c tan 0 . 


(34) 
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Similarly, by substituting from (23), (28), and Appendix A into (29) 
we obtain 

5 el 5e 

_r X- , jCOSy- T COS-y 

o 11 (x 1 ,x 2 ,x 3 ) = C^t 1 ) + X h 2 (1)] '^=7 ~ .2 yT * 

r s ° + -?-J 

Also if we express the (in-plane cleavage) stress in the close 
neighborhood of the crack tip (xi=0,x 2 =a) or (p=0) in the following 
standard form 


( 35 ) 


_ k,(x 3 ) 

a n ( x i ? x 2 ,x 3) = ^ f n (ct) ’ 


(36) 


with f i -j (0) =1, from (35) and (36) the Mode I stress intensity factor k-j 
and the function f-j-j giving the angular distribution of the stress 
component a-|-j around the crack tip may be obtained as follows 

x- 


k,(x 3 ) ■- i r[Vl) , T h l! | l) ] ' 


(37) 


5 0 1 50 

4 cos o ” 4 cos 2 1 

f n (a) = — — ; 2 - , e = Arctan(-tana) 


/„-2 L sin a\ 
^cos a + — 2 — J 


(38) 


Similarly if we let 


k-i (x 3 ) _ k, (x 3 ) 

a 22^ x l* x 2* x 3^ = f 22^* CT i2^ x r x 2 ,X 3^ ^ 


from (29) we obtain 


3 0^1 50 1 . 0 J 1 . 50 

f , . 4 cos^cos^r ._ sin _ +T sin T 

'22 a " 9 * » *19V a ) " 


/ 2 .sin^aY 1 * 

r 


12 


[ n .2 . sin 2 a^ 

(cos a+ ~J 


(39) 


(40) 


4. The Results and Discussion 

The symmetric problem formulated by the system of integral equations 
(23) is solved under two sets of external loads, namely uniform membrane 
loading and uniform bending moment applied to crack surfaces. For the 
membrane loading it is assumed that 
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N-j ] ( 0 * ) = ~ Nq = - h<? m » Wj ^ ( 0 ^ 2 ) = 0 , 


-a < x 2 < a , 


(41) 


where a m is constant. From (41) and the Appendix A the input functions 
in (23) may be expressed as 


F 1 ( y ) = " a n/ cE * F 2 M = 0 » -Sc<y<'ft- 


(42) 


In this problem the corresponding flat plate stress intensity factor is 
kp = c m v^r and from (37) the membrane and bending components of the stres 
intensity factor ratio k-j (x^J/o ^T may be evaluated as follows: 

k. 


mm 


= M 0) _ cE h m 

" " 2or m h l 1 * 

m 


m 


‘bm 


k^h/2) -^(0) cE h 


ay a 

m 


m 


For the bending of the shell the external loads are 

h 2 

N„(0,x 2 ) = 0 > H n (0,x 2 ) =-H 0 =- 5 » b . -a < x 2 < a 

which give the input functions as follows: 


(43) 


(44) 


(45) 


i r 1 (y) = o, F 2 (y) = ‘6 Ee» 


(46) 


where a b is a known constant. For this loading condition the membrane 
and bending components of the stress intensity factor ratio may be 
defined as and evaluated from 


k mb = 


M 0) _ cE hm 


bb 


k^h/2) - k, (0) E h 

‘“2^ 2iT h 2 (1) 




(47) 

(48) 


The calculated results for isotropic and specially orthotropic 
shells are shown in Figures 2-9 and Tables 1-11. Fig. 2 shows the com- 
parison of the membrane component of the stress intensity factor ratios 
k(nm ln a cylindrical shell containing a circumferential or an axial crack 
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and a spherical shell containing a meridional crack subjected to uniform 
membrane stress For the same loading and curvature the Mode I 

stress intensity factor appears to be highest in sphere and lowest for 
the circumferential crack. Similar results may be observed from Fig. 3 
which shows an example for the membrane stress resultant N-|j(0,x 2 ) for 
x 2 >a . After determining (*j and G 2 N^ is obtained directly from (23a) 
by observing that (23a) gives the expression of N xx (0,y) outside (i.e., 
for |y| >Jc) as well as inside the crack. Fig. 3 also shows the flat 
plate solution. Note that even though near the crack tip the stresses 
in the shells are greater than that in the plate, because of their greater 
rate of decay, away from the crack region the shell stresses fall below 
the stress level in the flat plate. 

The results showing the effect of the curvature ratio R|/R 2 on the 
Mode I stress intensity factor ratios k. • , (i,j=m,b) defined by (43), 

I J 

(44), (47) and (48) in an isotropic shell are given in Tables 1 and 2. 

Table 1 shows the results for a positive curvature ratio R-j/R 2 as, for 
example, in the case of outside surface of pipe elbows and barrel shaped 
toroidal shells. The results of an example for a negative curvature 
ratio R i / R 2 = 5 are given in Table 2. Fig. 4 shows the comparison of 

the stress intensity factors obtained from shells with positive and negative 
curvature ratios and from an axially cracked cylinder. Note that, as one 
might expect and consistent with the trends in Fig. 2, the stress intensity 
factor for R-j/R 2 >0 is greater and that for (R^/R 2 ) <0 is smaller than the 
value for the cylinder (i.e., for R-j/R^O). 

The remaining results given in this paper deal with the specially 
orthotropic shells, that is, the orthotropic shells for which the measured 
in-plane elastic constants of the material approximately satisfy the 
factorization condition given by (2). namely 



& 


1 E 2 


2(l+ify| v 2 ) 



(49) 


For example, consider two graphite-epoxy (fiber reinforced) laminates 
consisting of 0°, +45°, 90° unidirectional laminae. Laminate A has 20%, 
30%, and 50% of its laminae, and laminate B 20%, 50%, and 30% of its 
laminae oriented along 0°, +45°, and 90° directions, respectively. If 
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0° direction coincides with E-j , the measured elastic constants of the 


two laminates are 
E-j(psi) 

known to be 
E 2 (psi) 

V 1 

v 2 

G 12 (psi) 

G av (psi) 

A: 6.9 xio 6 

12.3 xio 6 

0.140 

0.250 

2.1 xio 6 

3.88 xio 1 

B: 7.1 xlO 6 

9 xio 6 

0.270 

0.342 

3.05 xio 6 

3.06 xio 1 


where 6 ay is calculated from (49) by using the effective modulus and the 
effective Poisson's ratio. Thus, it is seen that the assumption of special 
orthotropy G -|2 = G av is valid for laminate B but not for laminate A. 

Table 3 shows the elastic constants of two orthotropic materials used in 
the examples. The "mildly" orthotropic material is typical of rolled 
sheet metallic materials. Usually strongly orthotropic structural materials 
are fiber reinforced composites. 

Tables 4-9 show the membrane and bending components of the stress 
intensity factor ratio k mm and k bm defined by (43) and (44) for three 
symmetric crack geometries in cylindrical and spherical shells subjected 
to uniform membrane loading = - N on the crack surfaces (Fig. 1). 

For completeness the results for isotropic cylinders are also given in 
the tables. Note that for these simple crack-shell geometries there are 
three length parameters, namely mean radius R, thickness h, and half crack 
length a. Therefore, the solution must contain two dimensionless parameters 
which in these examples are assumed to be a/h and or x ^ which contains 
a//Rh (see Appendix A). In the tables the shell is designated by 
and the material is oriented in such a way that the crack is parallel to 
E 2 axis. The two sets of orthotropic results given im the same table 
correspond to a 90-degree material rotation. E-J/E 2 = 1 corresponds to the 
isotropic shell. These tables show that the influence of the material 
orthotropy on the stress intensity factors could be quite significant. 

A graphic demonstration of this effect in a cylindrical shell with an 
axial crack is shown in Fig. 5^. It is seen that the membrane component 

^ In Fig. 5 the shell parameter Xj^ = [l2(l-v 1 \) 2 ) ¥* a/v® is used as a 

basis in comparing kmm for different shells. For an accurate comparison 
one should have /vjvj = v(iso tropic) if is used, otherwise one should 
use a/v^Rh as a basis. The difficulty in such comparisons is compounded 
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of the stress intensity factor increases with increasing E-j/E 2 . This 
effect is quite definite and, for higher values of E-|/E 2 , 1S highly 
pronounced in cylindrical shells with an axial crack and in spherical 
shells. On the other hand in circumferentially cracked cylindrical shells 
the ortho tropy effect seems to be rather insignificant. It should be 
pointed out that the effect of material orthotropy such as that shown in 
Fig. 5 is not confined to shells. For example, similar results are ob- 
served in the plane elasticity solution of an orthotropic infinite strip 
containing a crack parallel to its boundaries [13]. In this problem if the 
stress intensity factor is plotted against a/H (2a and 2H being the crack 
length and strip width) with E-j/E 2 as the parameter, the result would be 
identical (in form) to that shown in Fig. 5, that is (in the terminology of 
Fig. 5) for Ej/E 2 >1 k-| would be greater and for E-|/E 2 <1 k-j would be 
smaller than the stress intensity factor for the corresponding isotropic 
plate. It should, of course, be noted that the stress intensity factor is 
a measure of the crack driving force in fracture analysis of the structural 
component. In comparing two materials (or two material orientations with 
respect to the crack plane) one must also consider the fracture resistance 
of the material. For example, in fiber reinforced composites having a 
crack in E 2 plane for E^/E 2 >1 not only the stress intensity factor but 
also the fracture resistance of the material would be expected to be 
greater than the corresponding values for a 90° rotated material (i.e., 
for E-j /E 2 < 1 ) . Hence from the analysis alone it is not possible to infer 
the fracture strength of orthotropic plates and shells. 

In a mildly orthotropic shell the effect of material orthotropy on 
the distribution of the membrane resultant N^-j(0,x 2 ) outside the crack is 
shown in Fig. 6. Near the crack tip the results shown in Fig. 6 are con- 
sistent with that of Fig. 5. However, as in Fig. 3, away from the crack 
tip the trend showing the effect of E^Eg on the stress distribution is 

also by the fact that /v^Vo or v enters Into the shell analysis 
independently as well as through X^. However, it has been shown that 
the effect of v on is rather insignificant [4]. Also, for the 
three materials used in Fig. 5, namely for (E^/E 2 ) =1, 1.38, 26.67 
the coefficient [l2(l-v^V2) ]** of a/ /Eh in X^ is 1.82, 1.84, 1.85, 
respectively. This means that in Fig. 5 even though the comparison is 
not made for shells with identical geometries, the error should not be 
very high and particularly the trend should be accurate. 
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reversed. 


From (36) -(40) it may be seen that not only the stress intensity 
factor but also the angular distribution of stresses around the crack 
tip is influenced by the material arthotropy. Fig. 7 shows an example 
for this effect in a strongly orthotropic material. Note that in isotropic 
materials the functions f..(a), (i,j=l,2) (see (36) and (39)) giving the 
angular distribution of stresses under Mode I loading condition are 
independent of the elastic constants (and are obtained from (36) and 
(39) by letting c=l). In the specially orthotropic materials f . . seems 

U. ■ J 

to depend on the stiffness ratio c = (E 1 /E 2 ) . 

Tables 4-9 show the stress intensity factors for shells under 
membrane loading only; i.e., for N^(0,x 2 ) = -N Q , M n (°>x 2 ) = 0 . For 
orthotropic shells subjected to bending on the crack surfaces (i.e., for 
N 11 (0,x 2 )*0, M-j-j (0,x 2 ) = -M q ) , the stress intensity factor ratios k bb and 
k mb de ‘ fr ' ined b y (47) and (48) do not seem to differ significantly from the 
corresponding results for the isotropic shells which are given in [3]-[5]. 
The results are, therefore, not tabulated in this paper. Some sample 
results are, however, shown in Figures 8 and 9 for a cylindrical shell 
containing an axial or a circumferential crack. The results for spherical 
shells may be found in [5]. 

Table 10 shows the principal stress intensity factor ratio k bb for 
a plate under bending. In addition to the relative crack size a/h, the 
table shows the effect of the Poisson's ratio in isotropic plates and 
the stiffness ratio E-j /E 2 in specially ortho tropic plates. Since the 
relative influence of the material constants on k bb in plates and shells 
is roughly the same, these results may be used with the isotropic shell 
results given for a fixed Poisson's ratio (e.g. , v=0.3 in most cases) to 
estimate the values of the k bb for the isotropic and orthotropic shells. 

The effect of material orthotropy on the stress intensity factor 
ratios k. ., (i,j=m,b) defined by (43), (44), (47) and (48) on a shell with 

* J 

principal radii of curvature R-j and R 2 is shown in Table 11 which also 
includes the isotropic results for comparison. Note that the results given 
in Table 11 have the same trend as that shown in Fig. 5 with regard to the 
variation of the principal stress intensity component k^ for the shell 
under membrane loading. 
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Appendix A 


The dimensionless and normalized quantities used in the analysis 


x = x-j/a/c , y = Xgv'c/a , z = x 3 /a ; (A.l) 

u = u-j/c/a , v = Ug/a^E, w = Ug/~a ; * (A. 2) 

3 x =3i^, B y = B 2 /^, <t»(x,y) = F(x 1 ,x 2 )/Eha 2 ; (A. 3) 

= * <7yy = G22 * a X y ^ °1 2/E » (A* 4) 

N xx = N n/ Ehc » N yy = cN 22 /Eh , N xy = N 12 /Eh ; (A. 5) 

M xx = M ll /Ech2 ’ M yy = cM 22 /Eh2 ’ M xy = M 12 /Eh2; (A - 6) 

V x = V-|/hB/c , V y = V 2 ^/Bh ; (A. 7) 

X 4 = 12(l-v 2 )a 4 c 2 /h 2 R 2 , X 4 = 12(l-v 2 )a 4 /c 2 h 2 R 2 , 

X 4 2 = 12(l-v 2 )a 4 /h 2 R 2 2 , X 4 = 12(l-v 2 )a 2 /h 2 , k = E/BX 4 ; (A.8) 

E = /E^\ v = (v ] /E 1 = v 2 /E 2 ) , B = 5E/1 2(1 +v) . (A.9) 
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Table 1. Stress intensity factor ratios in an isotropic shell with 

double curvature under uniform membrane load Nj-| =N 0 , Myj »0, 
( k mm» k bm) or uniform bending moment = M q , N n =0, U bb ,k mb ) 

(the crack in the 1^ plane, v = 0.3, a/h = 10) 


R 1 

a/ZhRg 

*2 

0.1 

0.25 

0.50 

0.75 

1.0 

1.5 

2.0 

k 

mm 

0.2 

IBSl 

1.216 

1.681 

2.246 

2.854 

4.145 

5.556 

1/3 

Wm 

1.141 

1.471 

1.897 

2.378 

3.434 

4.590 

0.5 

KBsm 

1.100 

1.349 

1.686 

2.078 

2.969 

3.965 

2 

1.005 

1.035 

1.130 

1.277 

1.460 

1.917 

2.471 

3 

1.004 

1.027 

1.101 

1.219 

1.368 

1.740 

2.195 

5 

1.003 

1.021 

1.078 

1.168 

1.286 

1.577 

1.934 




k bm 





0.2 

0.046 

0.134 

0.199 

0.142 

-0.038 

-0.742 

-1.899 

1/3 

0.034 

0.109 

0.188 

0.181 

0.081 

-0.392 

-1.221 

0.5 

0.028 

0.092 

0.174 

0.189 

0.132 

-0.207 


2 

0.015 

0.055 

0.123 

0.163 

0.166 

0.062 

HEbM] 

3 

0.013 

0.050 

0.112 

0.153 

0.161 

0.084 

fijiHIAj 

5 

0.011 

0.045 

0.102 

0.142 

0.153 

0.097 


k bb 

0.2 


0.615 

0.555 

0.495 


0.356 

0.293 

1/3 


0.621 

0.566 

0.508 


0.368 

0.305 

0.5 

0.644 

0.624 

0.573 

0.516 

0.463 

0.377 

0.313 

2 

0.645 

0.630 

0.586 

0.532 

0.480 

0.393 

0.331 

3 

0.645 

0.630 

0.587 

0.535 

0.482 

0.395 

0.334 

5 

0.645 

0.631 

0.589 

0.536 

0.484 

0.397 

0.336 

k mb 

0.2 

0.011 

0.033 

0.064 

0.086 

0.100 

0.115 

0.122 

1/3 

0.008 

0.026 

0.054 

0.074 

0.089 

0.104 

0.111 

0.5 

0.006 

0.022 

0.047 

0.066 

0.079 

0.095 


2 

0.003 

0.013 

0.030 

0.044 

0.055 

0.067 

0.073 

3 

0.003 

0.011 

0.027 

0.040 

0.050 

0.061 


5 

0.003 

0.010 

0.024 

0.037 

0.045 

0.054 
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Table 2. Stress intensity factors in a saddle-shaped shell under uniform 
membrane loading N„=N 0 , Mp-O, ^ k mm’ k bm^‘ or uniform bending 
M^I =M 0> N-j^ =0, (the crack in plane, R^/R 2 = -0.5, 

a/h = 10, v = 0.3) 


a/ZRgh 

0.1 

— 

0.25 

0.5 

0.75 

1.0 

1.5 

2.0 

Snm 

1.014 

1.079 

1.261 

1.480 

1.699 

2.098 

2.455 

k bm 

-0.015 

-0.045 

-0.066 

-0.039 

0.025 

0.195 

0.335 

k bb 

0.645 

0.633 

0.594 

0.546 

0.498 

0.414 

0.353 

k mb 

-0.003 

-0.011 

-0.019 

-0.020 

-0.016 

-0.004 

0.006 

1 


Table 3. Elastic constants of the orthotropic materials used in the 
examples (v-j/E-j = v 2 /E 2> G av = E/2(l+v) , E = /E^E^", v = 



Mildly ortho tropic 
material (titanium) 

Strongly orthotropic 
material (graphite-epoxy) 

E 1 N/m 2 (psi) 

11 7 

1.039 *10" (1 .507 * 10 ; 

1.034 xlO 10 (1.5 xio 6 ) 

E 2 N/m 2 (psi) 

1.434 xlO 11 (2.08 xlO 7 ) 

2.758 xlO 11 (4 xlO 7 ) 

V 1 

0.1966 

0.0075 

v 2 

0.2714 

0.2000 

G 12 N/m 2 (psi) 

4.675 xlO 10 (6.78 xlO 6 ) 

2.758 xlO 10 (4.0 xlO 6 ) 

G av N / m2 ( psi ) 

4.930 xlO 10 (7.15 xlO 6 ) 

2.572 xlO 10 (3.73 xlO 6 ) 
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Table 4. Stress intensity factors in a mildly orthotropic and in an isotropic cylindrical shell containing 
an axial through crack and subjected to uniform membrane loading N 11= N 0 (v = 0.3 for the isotropic 
shell x 1 = [12(1 -v 2 )] 1 * a/ M, v = /y£) 



10 

5 

2 

0.725 1.380 1 

0.725 1.380 1 

0.725 1.380 1 



0.725 1.380 

















Table 5. Stress intensity factors in a strongly orthotropic and in an isotropic (v = 0.3) cylindrical 
shell containing an axial crack and subjected to uniform membrane loading =N Q . 


a/h 

1° 

5 

2 

1 

e,/e 2 

0.037 26.67 

1 

0.037 26.67 

1 

0.037 26.67 

1 

0.037 26.67 

1 


X 1 k mm 


0.0 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

0.5 

1.012 

1.238 

1.057 

1.012 

1.238 

1.057 

1.012 

1.242 

1.058 

1.012 

1.253 

1.061 

1.0 

1.043 

1.714 

1.199 

1.043 

1.717 

1.200 

1.045 

1.736 

1.208 

1.048 

* 

1.795 

1.233 

1.5 

1.093 

2.244 

1.394 

1.093 

2.252 

1.398 

1.047 

2.298 

1.420 



2.0 

1.155 

2.779 

1.618 

1.156 

2.794 

1.625 

1.164 

2.877 

1.668 




3.0 

1.309 

3.826 

2.105 

1.313 

3.865 

2.122 







4.0 

1.490 

4.852 

2.603 

1.498 

4.955 

2.634 







5.0 

1.686 

— 

3.096 

1.699 

— 

3.146 







6.0 

1.891 

— 

3.580 










8.0 

2.311 

— 

4.515 







1 



10.0 

2.737 

— 

5.422 











k bm 


0.0 



0.000 


lj|K|i|i| 


0.000 

nn 


0.000 


ill 

0.5 

1 ffiVM 

Kfjjt m 

0.054 

1 



0.009 

^ Mi 


0.010 



1.0 


■nftj m 

0.120 


■IWtH 


0.024 


mmnm 

0.027 


warn 

1.5 

■S)(0cf*fl 


0.167 

0.040 

0.172 

0.164 

0.042 

0.151 

0.160 




2.0 

0.056 

0.175 

0.194 

0.056 

0.148 

0.185 

0.059 

0.105 

0.175 




3.0 

0.089 

0.017 

0.189 

0.088 

-0.045 

0.166 







4.0 

0.117 

-0.287 

0.116 

0.115 

-0.395 

0.076 







5.0 

0.141 

— 

-0.014 

0.136 

— 

-0.070 







6.0 

0.158 

— 

-0.187 










8.0 

0.169 

— 

-0.634 










10.0 

0.148 

— 

-1.183 












Table 6. Stress Intensity factors in a mildly orthotropic and in an isotropic cylindrical shell 
containing a circumferential through crack and subjected to uniform membrane loading 
Nn =N 0 , ( v = 1/3 for the isotropic shell, a 2 = [12(l-v 2 )] Jzi a/M\, v = /v-j v 2 ) 


a/h 

10 

5 

2 

1 

e,/e 2 

0.725 1.380 

1 

0.725 1.380 

1 

0.725 1.380 

1 

0.725 1.380 

1 


X 2 k mm 


0.0 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

0.5 

1.012 

1.012 

1.012 

1.012 

1.012 

1.012 

1.013 

1.013 

1.012 

1.013 

1.013 

1.013 

1.0 

1.047 

1.047 

1.048 

1.047 

1.047 

1.048 

1.049 

1.049 

1.050 

1.053 

1.053 

1.055 

1.5 

1.100 

1.099 

1.102 

1.101 

1.100 

1.103 

1.105 

1.104 

1.108 




2.0 

1.164 

1.164 

1.168 

1.166 

1.165 

1.169 

1.174 

1.173 

1.179 




3.0 

1.308 

1.308 

1.314 

1.312 

1.311 

1.317 








1.457 

1.456 

1.462 

1.462 

1.461 

1.467 








1.599 

1.598 

1.604 











k bm 



0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


0.000 

0.000 

nn 1 

0.036 

0.036 

0.041 

0.036 

0.036 

0.041 

0.038 

0.038 

0.043 

1 MB 

0.040 

0.044 


0.081 

0.082 

0.092 

0.081 

0.081 

0.092 

0.082 

0.082 

0.092 

■tlliMW 

0.083 

0.093 

1.5 

0.110 

0.111 

0.123 

0.108 

0.109 

0.119 

0.104 

0.105 

0.114 





0.115 

0.117 

0.125 

0.110 

0.111 

0.119 

0.100 

0.101 

0.107 





0.072 

0.075 

0.071 

0.059 

0.062 

0.057 








-0.009 

-0.006 

-0.024 

-0.026 

-0.024 

-0.042 








-0.100 

-0.096 

-0.126 

















Table 7. Stress intensity factors in a strongly orthotropic and in an isotropic ( v = 1/3) cylindrical 
shell containing a circumferential crack and subjected to uniform membrane loading =N Q . 


a/h 

10 

5 

2 

1 

m 

m 

ro 

0.037 26.67 1 

0.037 26.67 1 

0.037 26.67 1 

0.037 26.67 1 



0.0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.5 

0.028 

0.028 

0.041 

0.029 

0.028 

0.041 

0.030 

0.029 

0.043 

0.033 

0.030 

0.044 

1.0 

0.066 

0.065 

0.092 

0.066 

0.065 

0.092 

0.067 

0.065 

0.092 

0.069 

0.066 

0.093 

1.5 

0.090 

0.091 

0.123 

0.088 

0.089 

0.119 

0.084 

0.087 

0.114 




2.0 

0.096 

0.101 

0.125 

0.090 

0.099 

0.119 

0.080 

0.090 

0.107 




3.0 

0.065 

0.077 

0.071 

0.051 

0.067 

0.057 







4.0 

0.001 

0.019 

-0.024 

-0.017 

0.004 

-0.042 







5.0 

-0.072 

-0.050 

-0.126 












Table 8. Stress intensity factors in a mildly orthotropic and in an isotropic (v = l/3) spherical shell 
containing a meridional crack and subjected to uniform membrane loading = N Q 
(x £ = [12(1 -v 2 )]^ a//Rh, v = 


a/h 

10 

5 

2 

1 

V E 2 

0.725 1.380 

1 

0.725 1.380 

1 

0.725 1.380 

1 | 

| 0.725 1.380 

1 


x 0 k 

2 mm 


0.0 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1 .000 

1.000 

1.000 

1.000 

1.000 

1.000 

0.25 

1.016 

1.021 

1.018 

1.016 

1.021 

1.018 

1.016 

1.021 

1.019 , 

1.017 

1.022 

1.019 

0.50 

1.060 

1.077 

1.069 

1.060 

1.077 

1.069 

1.062 

1.079 

1.071 ! 

1.065 

1.084 

1.076 

0.75 

1.129 

1.164 

1.149 

1.130 

1.165 

1.150 

1.135 

1.171 

1.156 

1.148 

1.187 

1.173 

1.0 

1.220 

1.276 

1.252 

1.222 

1.279 

1 .255 

1 .232 

1.291 

1.268 

1.262 

1.327 

1.305 

1.5 

1.450 

1.555 

1.512 

1.456 

1.561 

1.519 

1.486 

1.596 

1.556 




2.0 

1.733 

1.891 

1.828 

1.744 

1.903 

1.841 

1.807 

1.976 

1.918 | 




2.5 

2.056 

2.270 

2.186 

2.076 

2.291 

2.208 







3.0 

2.415 

2.685 

2.579 

2.446 

2.718 

2.615 



l 




3.5 

2.806 

3.132 

3.004 

2.852 

3.183 

3.058 







4.0 

3.226 

3.611 

3.460 

3.296 

3.687 

3.539 








k bm 


0.0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.25 

0.025 

0.029 

0.033 

0.026 

0.030 

0.034 

0.027 

0.031 

0.035 

0.029 

0.033 

0.037 

0.50 

0.065 

0.073 

0.084 

0.066 

0.074 

0.084 

0.068 

0.075 

0.086 

0.071 

0.079 

0.090 

0.75 

0.104 

0.115 

0.130 

0.104 

0.115 

0.130 

0.106 

0.115 

0.130 

0.109 

0.118 

0.133 

1.0 

0.135 

0.148 

0.165 

0.134 

0.145 

0.162 

0.133 

0.142 

0.158 

0.134 

0.141 

0.157 

1.5 

0.165 

0.175 

0.187 

0.158 

0.164 

0.174 

0.145 

0.146 

0.155 




2.0 

0.146 

0.145 

0.142 

0.129 

0.121 

0.117 

0.099 

0.083 

0.077 




2.5 

0.077 

0.054 

0.031 

0.047 

0.015 

-0.010 







3.0 

-0.044 

-0.099 

-0.146 

-0.090 

-0.158 

-0.206 







3.5 

-0.219 

-0.317 

-0.390 

-0.283 

-0.398 

-0.471 







4.0 

-0.448 

-0.602 

-0.701 

-0.533 

-0.709 

-0.807 









Table 9. Stress intensity factors in a strongly orthotropic and in an isotropic (v = 1/3) spherical 
shell containing a meridional crack and subjected to uniform membrane loading = N Q . 


a/h 

10 

5 

2 

1 

m 

m 

ro 

0.037 26.67 

1 

0.037 26.67 

1 

0.037 26.67 

1 

0.037 26.67 

1 


0.0 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

0.25 

1.006 

1.073 

1.018 

1.006 

1.073 

1.018 

1.006 

1.074 

1.019 

1.006 

1.075 

1.019 

0.50 

1.024 

1.250 

1.069 

1.024 

1.251 

1.069 

1.024 

1.255 

1.071 

1.026 

1.267 

1.076 

0.75 

1.051 

1.491 

1.149 

1.052 

1.493 

1.150 

1.053 

1.504 

1.156 

1.059 

1.540 

1.173 

1.0 

1.089 

1.767 

1.252 

1.090 

1.771 

1.255 

1.095 

1.794 

1.268 

1.108 

1.863 

1.305 

1.5 

1.190 

2.377 

1.512 

1.193 

2.387 

1.519 

1.207 

2.445 

1.556 




2.0 

1.318 

3.030 

1.828 

1.324 

3.050 

1.841 

1.354 

3.165 

1.918 




2.5 

1.470 

3.713 

2.186 

1.480 

3.748 

2.208 







3.0 

1.643 

4.423 

2.579 

1.658 

4.484 

2.615 







3.5 

1.835 

5.164 

3.004 

1.857 

5.271 

3.058 







4.0 

2.044 

5.948 

3.460 

2.075 

6.145 

3.539 








k bra 


0.0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


BH 


{Vi fflj 

liWiTWl 

0.000 

0.25 

0.012 

0.037 

0.033 

0.012 

0.037 

0.034 




1 


0.037 

0.50 

0.033 

0.088 

0.084 

0.033 

0.087 

0.084 






0.090 

0.75 

0.055 

0.133 

0.130 

0.056 

0.130 

0.130 


0.126 

0.130 

Wmi inB 

MmnM 

0.133 

1.0 

0.076 

0.165 

0.165 

0.076 

0.158 

0.162 


0.147 

0.158 

0.079 

0.141 

0.157 

1.5 

0.104 

0.172 

0.187 

0.101 

0.151 

0.174 


0.114 

0.155 




2.0 

0.111 

0.084 

0.142 

0.103 

0.041 

0.117 

*5 

-0.035 

0.077 




2.5 

0.099 

-0.108 

0.031 

0.085 

-0.179 

-0.010 







3.0 

0.069 

-0.404 

-0.146 

0.048 

-0.510 

-0.206 







3.5 

0.022 

-0.804 

-0.390 

-0.006 

-0.957 

-0.471 







4.0 

-0.041 

-1.314 

-0.701 

-0.074 

-1.538 

-0.807 










Table 10. The stress intensity factor ratio in isotropic and 

specially orthotropic plates under bending (the crack along 
the Xg axis, v = /v-j v-j , Ej/E 2 = 1.38, 0.725 for titanium, 
E-j/Eg 3 26.67, 0.037 for graphite epoxy) 


h'h 

V 

a/h = 10 

a/h = 5 

a/h = 2 

a/h = 1 

1 

0 

0.598 

0.613 

0.651 

0.702 

1 

0.1 

0.616 

0.631 

0.669 

0.719 

1 

0.2 

0.632 

0.647 

0.685 

0.734 

1 

0.3 

0.644 

0.662 

0.699 

0.747 

1 

1/3 

0.652 

0.667 

0.704 

0.752 

1 

0.4 

0.661 

0.675 

0.712 

0.760 

1 

0.5 

0.673 

0.688 

0.724 

0.771 

0.725 

0.231 

0.632 

0.651 

0.691 

0.742 

1.38 

0.231 

0.634 

0.650 

0.687 

0.735 

0.037 

0.037 

0.615 

0.635 

0.686 

0.751 

26.67 

0.037 

0.599 

0.611 

0.639 

0.677 
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Table 11. Stress intensity factor ratios in an isotropic and in a specially orthotropic shell with 

principal radii of curvature Rj and Rg having a through crack in R£ plane (material: titanium, 
R-|/R2 = l/3, a/h = l/10, for the isotropic shell v = 0.3) 




a/Z^F 


h'h 

0.1 

0.25 

0.50 

0.75 

1.0 

1.5 

2.0 


0.725 

1.022 

1.122 

1.412 

1.795 

2.233 

3.210 

4.294 

k mm 

1.38 

1.029 

1.159 

1.521 

1.981 

2.492 

3.604 

4.824 

1.0 

1.025 

1.141 

1.471 

1.897 

2.378 
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Fig. 2 The comparison of the membrane components of the stress intensity 
ratio = for a cylindrical shell with a circumferential 

or an axial crack and a spherical shell with a meridional crack 
all under uniform membrane stress a . 
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Fig. 3 Distribution of the membrane resultant li^ outside the crack in 
isotropic spherical and cylindrical shells and flat plates 
v = 0.3, a/h =5, X = [12(1 -v 2 ) a//Rh = 1.5 in shells and X = 0 
in plates. 
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Fig. 5 The influence of material orthotropy on the membrane component 
of the stress intensity factor in a cylindrical shell containing 
an axial crack; (a/h = 10, v = 0.3 for the isotropic shell) 



Fig. 6 The influence of the material orthotropy on the distribution of 
membrane resultant 1^(0, Xg) outside the crack in a cylindrical 
shell containing an axial crack; (a/h = 5, A-j =1.5, v = 0.3 for 
the isotropic shell) 
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Fig. 7 The influence of the material orthotropy on the angular distribu- 
tion 1 1 ( “ ) of the stress component for Mode I loading 

condition (the crack lies along the axis) 
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Fig. 8 The influence of the material orthotropy on the membrane and 

bending components of the stress intensity factor in a cylindrical 
shell containing an axial crack and subjected to uniform bending 
(0,X£) = -M q on the crack surfaces; (a/h = 5, v = 0.3 for the 
isotropic shell, x = [12(1 -v 2 ) a//Rh) 
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Fig. 9 Same as Fig. 8 for a cylindrical shell with a ci 
crack; (a/h = 5, v = l/3 for the isotropic shell) 
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